A Signature of Cosmic Strings Wakes in the CMB Polarization 
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We calculate a signature of cosmic strings in the polarization of the cosmic microwave background 
(CMB). We find that ionization in the wakes behind moving strings gives rise to extra polarization in 
a set of rectangular patches in the sky whose length distribution is scale-invariant. The length of an 
individual patch is set by the co-moving Hubble radius at the time the string is perturbing the CMB. 
The polarization signal is largest for string wakes produced at the earliest post-recombination time, 
and for an alignment in which the photons cross the wake close to the time the wake is created. 
The maximal amplitude of the polarization relative to the temperature quadrupole is set by the 
overdensity of free electrons inside a wake which depends on the ionization fraction / inside the 
wake. The signal can be as high as 0.06/iK in degree scale polarization for a string at high redshift 
(near recombination) and a string tension ^ given by = 10~^. 
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PACS numbers: 98.80.Cq 



I. INTRODUCTION 



In recent years there has been renewed interest in the 
possibility that cosmic strings contribute to the power 
spectrum of curvature fluctuations which give rise to the 
large scale structure and cosmic microwave background 
(CMB) anisotropics which we see today. One of the rea- 
sons is that many inflationary models constructed in the 
context of supergravity models lead to the formation of 
gauge theory cosmic strings at the end of the inflation- 
ary phase [v\ . Secondly, in a large class of brane inflation 
models the formation of cosmic superstrings at the 
end of inflation is generic 0], and in some cases (see [1]) 
these strings are stable (see also [s'l for reviews on fun- 
damental cosmic strings). Cosmic superstrings are also 
a possible remnant of an early Hagedorn phase of string 
gas cosmology Q 

In models which admit stable strings or superstrings, 
a scaling solution of such strings inevitably results as 
a consequence of cosmological dynamics (see e.g. [1] for 
reviews on cosmic strings and structure formation). In 
a scaling solution, the network of strings looks statisti- 
cally the same at any time t if lengths are scaled to the 
Hubble radius at that time. The distribution of strings is 
dominated by a network of "inflnite" strings ^ with mean 
curvature radius and separation being of the order of the 
Hubble radius. The scaling solution of infinite strings is 
maintained by the production of string loops due to the 
interaction of long strings. This leads to a distribution 
of string loops with a well defined spectrum (see e.g. 0) 
for all radii R smaller than a cutoff radius set by the 
Hubble length. Whereas the scaling distribution of the 
infinite strings is reasonably well known as a result of 
detailed numerical simulations of cosmic string evolution 



^ Any string with a mean curvature radius comparable or greater 
than the Hubble radius or which extends beyond the Hubble 
radius is called "infinite" or "long" . 



(see [13 for some references), there is still substantial 
uncertainty concerning the distribution of string loops. 
It is, however, quite clear that the long strings dominate 
the energy density of strings. The distinctive signals of 
strings which we will focus on are due to the long strings. 

Cosmic strings give rise to distinctive signatures in 
both the CMB and in the large-scale structure. These 
signatures are a consequence of the specific geometry of 
space produced by strings. As studied initially in [1 1| . 
space perpendicular to a long straight string is locally 
flat but globally looks like a cone whose tip coincides 
with the location of the string (the smoothing out of the 
cone as a consequence of the internal structure of the cos- 
mic string was worked out in The deficit angle is 
given by 



a = SvrG/i . 



(1) 



where /i is the string tension and G is Newton's con- 
stant. Hence, a cosmic string moving with velocity v in 
the plane perpendicular to its tangent vector will lead to 
line discontinuities in the CMB temperature of photons 
passing on different sides of the string. The magnitude 
of the temperature jump is [ij] 



5T_ 



(2) 



where 7 is the relativistic gamma factor associated with 
the velocity v. 

As a consequence of the deficit angle H]), a moving 
string will generate a cosmic string wake, a wedge-shaped 
region behind the string (from the point of view of its 
velocity), a region with twice the ba ckg round density [14 1 
(see Figure 1). Causality (see e.g. [l^) limits the depth 
of the distortion of space due to a cosmic string. The 
details were worked out in [l^ where it was shown that 
the deficit angle goes to zero quite rapidly a distance t 
from the string. Hence, the depth of the string wake is 
given by the same length. String wakes lead to distinctive 
signatures in the topology of the large-scale structure, 
signatures which were explored e.g. in [l7j . 
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FIG. 1: Sketch of the mechanism by which a wake behind a 
moving string is generated. Consider a string perpendicular 
to the plane of the graph moving straight downward. From 
the point of view of the frame in which the string is at rest, 
matter is moving upwards, as indicated with the arrows in 
the left panel. From the point of view of an observer sitting 
behind the string (relative to the string motion) matter flow- 
ing past the string receives a velocity kick towards the plane 
determined by the direction of the string and the velocity vec- 
tor (right panel). This velocity kick towards the plane leads 
to a wedge-shaped region behind the string with twice the 
background density (the shaded region in the right panel). 



Wakes formed at arbitrarily early times are non-linear 
density perturbations. For wakes formed by strings 
present at times U > tree, where tree is the time of re- 
combination, the baryonic matter inside the wake un- 
dergoes shocks [3 (see e.g. [l^ for a detailed study). 
The shocks, in turn, can ionize the gas - although as 
it turns out the residual ionization from decoupling is 
larger. Photons passing through these ionized regions on 
their way from the last scattering surface to the observer 
can thus be polarized - and it is this polarization signa- 
ture which we aim to study here. 

The tightest constraints on the contribution of scal- 
ing strings to structure formation (and thus the tight- 
est upper bound on the tension fi of the strings) comes 
from the analysis of the angular power spectrum of 
CMB anisotropics. As discussed in |20l-|22|. the angular 
power spectrum does not have the acoustic ringing which 
inflation-seeded perturbations generate. The reason is 
that the string network is continuously seeding the grow- 
ing mode of the curvature fluctuation variable on super- 
Hubble scales. Hence, the fluctuations are "incoherent" 
and "active" as opposed to "coherent" and "passive" as 
in the case of inflation-generated fluctuations. The con- 
tribution of cosmic strings to the primordial power spec- 
trum of cosmological perturbations is thus bounded from 
above, thus leading to an upper bound [23l - l27j on the 
string tension of about G/i < 3 x 10^^. 

Past work on CMB temperature maps has shown that 
signatures of cosmic strings are easier to identify in po- 
sition space than in Fourier space [l^, [28l - [30[ and recent 
studies show that high angular resolution surveys such as 
the South Pole Telescope project [3l| have the potential 



of improving the limits on the string tension by an order 
of magnitude [32-34]. 

To date there has been little work on CMB polariza- 
tion due to strings. Most of the existing work focuses 
on the angular power spectra of the polarization. Based 
on a formalism (see also [s^) to include cosmic de- 
fects as source terms in the Boltzmann equations used 
in CMB codes, the power spectra of temperature and 
temperature polarization maps were worked out 36] in 
the case of models with global defects such as global cos- 
mic strings. Since in cosmic defect models vector and 
tensor modes are as important as the scalar metric fluc- 
tuations [2^ \3T\ , a signiflcant B-mode polarization is in- 
duced. In fact, in the case of global strings with a tension 
close to the upper bound mentioned above, whereas the 
amplitude of the temperature and E-mode polarization 
power spectra are so small as to make the string signal 
invisible compared to the signal from the scale-invariant 
spectrum of adiabatic fluctuations (e.g. produced in in- 
flation) , the contribution of strings dominates the ampli- 
tude of the B-mode polarization power spectrum. In the 
case of local strings, these conclusions were confirmed in 
the more recent analyses of H^, [2^, [l^, [H, . The max- 
imal amplitude of the B-mode polarization power spec- 
trum for strings with G/i = 3 x 10~^ was shown to be 
taken on at angular harmonic values of I ~ 500 and to 
be of the order O.S^iK^ [39|]. However, the analyses of 
[H m [13, m, Hi] do not take into account the effects 
of the gravitational accretion onto cosmic string wakes. 
In related work, the conversion of E-mode to B-mode po- 
larization via the gravitational lensing induced by cosmic 
strings was studied in ^. 

Similarly to what was found in the analysis of CMB 
temperature maps from cosmic strings, we expect that a 
position-space analysis will be more powerful at revealing 
the key non- Gaussian signatures of strings in CMB polar- 
ization. Hence, in this work we derive the position space 
signature of a cosmic string wake in CMB polarization 
maps. 

In the following section we briefly review cosmic string 
wakes. In Section 3 we then analyze the polarization sig- 
nature of wakes, and we conclude with some discussion. 



II. COSMIC STRING WAKES 

Since the strings are relativistic, they generally move 
with a velocity of the order of the speed of light. There 
will be frequent intersections of strings. The long strings 
will chop off loops, and this leads to the conclusion that 
the string distribution will be statistically independent 



^ While this paper was being finalized for submission, a preprint 
appeared [4111 computing the local B-mode polarization power 
spectrum from cosmic strings. 
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on time scales larger than the Hubble radius 

In this work, we place one string of length citi ^ at 
a specified time ti and assume it is moving in transverse 
direction with a velocity Vs ■ This string segment will gen- 
erate a wake, and it is the signal of one of these wakes in 
the CMB polarization which we will study in the follow- 
ing. 

A string segment laid down at time ti. will generate a 
wake whose dimensions at that time are the following: 



CiU X tiVs^s X ^TtG^UVs^s 



(3) 



In the above, the first dimension is the length in direction 
of the string, the second is the depth which depends on 
Vs (7s is the associated relativistic 7 factor) , and the 
third is the average thickness. 

Once the wake is formed, its planar dimensions will 
expand as the universe grows in size, and the thickness 
(defined as the region of non-linear density) will grow by 
gravitational accretion. The accretion of matter onto a 
cosmic string wake was studied in ji ^ [43 in the case of 
the dark matter being cold, and in 44? 45| in the case 
of the dark matter being hot ^ . We are interested in the 
case of cold dark matter. 

We consider mass planes at a fixed initial comoving dis- 
tance q above the center of the wake. The corresponding 
physical height is 



h{q,t) = a{t)[q-4,{q,t)\ , 



(4) 



where ip{q,t) is the comoving displacement induced by 
the gravitational accretion onto the wake. For cold dark 
matter, the initial conditions for ijj{q) are il){q,ti) = 
''Pi'l: ti) — 0- The goal of the analysis is to find the thick- 
ness of the wake at all times t > ti. The thickness is 
defined as the physical height h above the center of the 
wake of the matter shell which is beginning to fall towards 



^ Hence, to model the effects of strings we will make use of a toy 
model introduced in f42| and used in most analytical work on 
cosmic strings and structure formation since then: we divide the 
time interval between recombination and the current time to into 
Hubble expansion time steps. In each time step, there is a dis- 
tribution of straight string segments moving in randomly chosen 
directions with velocities chosen at random between and 1 (in 
units of the speed of light). The centers and directions of these 
string segments are random, and the string density corresponds 
to N strings per Hubble volume, where N is an integer which 
is of the order 1 according to the scaling of the string network. 
The distribution of string segments is uncorrelated at different 
Hubble times. 

The constant ci is of the order 1 and depends on the correlation 
length of the string network as a function of the Hubble radius 
and must be determined from numerical simulations of cosmic 
string evolution. 

® If the strings have lots of small-scale structure then they will 
have an effective tension which is less than the effective energy 
density [4^ . This will lead to a local gravitational attraction 
of matter towards the string, a smaller transverse velocity, and 
hence to string filaments instead of wakes. The gravitational 
accretion onto string filaments was studied in [47| . 



the wake, i.e. for which h{q,t) = 0. In the Zel'dovich 
approximation, we first consider the equation of motion 
for h obtained by treating the source (the initial surface 
density cr of the wake) in the Newtonian limit, i.e. 



^= -dh 



(5) 



where $ is the Newtonian gravitational potential given 
by the Poisson equation 



92$ 



4TTG[p + a5{h)] 



(6) 



{p{t) being the background energy density), and then by 
linearizing the resulting equation in ip. The mean surface 
density is 



ait) - 47TGpt,Vsjs{lY^^pit)^ 



(7) 



The result of the computation of the value of the co- 
moving displacement qni which is "turning around" at 
the time t for a wake laid down at time U is 1441 



with 



qni{t,ti) = ^0(7)^^^ 



247r 



(8) 



MU) = —GpVsis{z{ti) + l)-^'ho. (9) 
5 

This corresponds to a physical height of 

h{t,ti) = a{t)[qni{t,ti) -ip{qni,t)] ~ a{t)qni{t,ti) 

where Zi and z are the redshifts corresponding to the 
times ti and t, respectively. These formulas agree with 
what is expected from linear cosmological perturbation 
theory: the fractional density perturbation should in- 
crease linearly in the scale factor which means that the 
comoving width of the wake must grow linearly with a{t). 

Let us return to the geometry of the string segment. 
The tangent vector to the string and the direction of 
motion of the string determine a two dimensional hyper- 
surface in space (the "string plane"). If neither the string 
tangent vector nor the velocity vector have a radial com- 
ponent (radial with respect to the co-moving point of our 
observer), then the photons will cross each point of the 
string wake at the same time, and the wake will corre- 
spond to a rectangle in the sky whose planar dimensions 
are 



ciii X t.Vsjs 



(11) 



However, if the normal vector of the string plane is not 
radial, then one of the planar dimensions will be reduced 
by a trigonometric factor cos(0) depending on the angle 
9 between the normal vector and the radial vector. In 
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addition, photons which we detect today did not pass the 
wake at the same time. This wih lead to a gradient of the 
polarization signal across the projection of the wake onto 
the CMB sky. To simplify the analysis, in the following 
we will assume that 9 is small. 

In the following section we will need the expression for 
the number density ne{t^ti) of free electrons at time t in 
a wake which was laid down at the time ti. The initial 
number density is 



(12) 



where / is the ionization fraction, pB is the energy den- 
sity in baryons, and mp is the proton mass. Taking into 
account that for t > ti the number density redshifts as 
the inverse volume we get 



(13) 



III. ANALYSIS 

If unpolarized cosmic microwave background radiation 
with a quadrupolar anisotropy scatters off free electrons, 
then the scattered radiation is polarized (see e.g. for 
reviews of the theory of CMB polarization). The magni- 
tude of the polarization depends on the Thomson cross 
section gt and on the integral of the number density of 
electrons along the null geodesic of radiation. Our start- 
ing formula is 



P(n) 



-(-) 



1/2 



(14) 



where P(n) is the magnitude of polarization of radiation 
from the direction n, Q is the temperature quadrupole, 
and 



TT = <7T i ne{x)dx, 



(15) 



where the integral is along the null geodesic in terms of 
conformal time. 

Let us now estimate the contribution to the polariza- 
tion amplitude from CMB radiaton passing through a 
single wake at time t, assuming that the wake was laid 
down at the time U and has thus had time to grow from 
time ti to the time t when the photons are crossing it. 
We assume that the photons cross in perpendicular direc- 
tion. Since the wake is thin, we can estimate the integral 
in Eq. ^ by 



TT 



2aTne{t,t,){z{t) + l)h{t,t^) 



(16) 



where the factor 2 is due to the fact that the width of 
the wake is twice the hight, and the redshift factor comes 
from the Jacobean transformation between x and posi- 
tion. Inserting Eq. (|T6| into the expression p4|) for the 
polarization amplitude and using the value for the height 



and the number density of electrons obtained in the pre- 
vious section we find 



P 

Q 



li^y^^'^TfpB{ti)m/ 

If 



5 Mtt 
{zjt) 



r -00(^1 



(17) 



Inserting the formula for V'o(^i) from Eq. (jH]), expressing 
the baryon density at ti in terms of the current baryon 
density, and the latter in terms of the baryon fraction 
^Ib and the total energy density pc at the current time 
to, we obtain 



"aTfGpVs^is 



P ^ 247r . 3 ^^l/2 
Q ~ "25"'4^ 

xl7BPc(to)m;Ho(z(i) -f l)'(z(<0 + 1) 



1)'/'(18) 



From Equation we see that the polarization signal 
is larger for wakes laid down early, i.e. close to the time of 
recombination. For fixed ti, the signal is largest for con- 
figurations where the photons we observe today cross the 
wake at the earliest possible time, i.e. for the largest z{t) 
(obviously, t is constrained to be larger than ti, other- 
wise our formula for the height is not applicable). To get 
an order of magnitude estimate of the magnitude of the 
polarization signal, we take z{ti) ~ z(t) ~ 10'^. Inserting 
the value of the Thomson cross section, the proton mass 
and the current time we get 



103 



1 



103 



)^o^ (19) 



The ionization fraction of baryonic matter drops off af- 
ter recombination, but it does not go to zero. As already 
discussed in [H, [FO] there is remnant residual ionization 
of the matter. As computed in [505 (see also '5l'|), the 
residual ionization fraction / tends to a limiting value of 
between 10^^ and lO""' at late times after recombination. 

Shocks inside the wake will lead to extra ionization. 
However, the resulting contribution to the ionization 
fraction is negligible for the range of string tensions we 
are interested in. To see this, we follow the analysis of 
[l8| . A particle streaming towards a cosmic string wake 
has velocity 



Vi = inGpVsjs 



(20) 



and kinetic energy ^ravf . The particles will undergo 
shocks and thermalize. Equating the initial kinetic en- 
ergy density with the final thermal energy density (as- 
suming that the particles are distributed according to an 
approximate ideal gas law) yields a temperature inside 
the wakes of 



r ~ 7 X 103 



Gp 



2 X 10- 



[Vslsf K . 



(21) 



We then compute the Boltzmann factor to determine 
the ratio of ionized Hydrogen to ground state Hydro- 
gen which is the ionization fraction. For string tensions 
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of order 10^^ this results in ionization fractions of the 
order 10^^ which is considerably less than the residual 
ionization. For more realistic string tensions compatible 
with current bounds, 10~^, the ionization fraction due to 
shocks is negligible compared to the residual ionization 
fraction. At the temperatures considered for a string 
tension of 10~^ at lower redshifts there would be star 
formation and hence an ionization fraction due to that 
effect since the wake would contain molecular hydrogen 
and satisfy the appropriate conditions. However, this is 
not an issue at redshifts under consideration. 

Note that even though the extra ionization inside the 
wake is negligible compared to the overall ionization level, 
the wake is a locus of extra energy density. Thus, even if 
the ionization fraction is homogeneous in space, the in- 
homogeneous distribution of matter will lead to a specific 
polarization signature. 

The position space signal of the polarization produced 
by a wake will be very specific; for a wake whose nor- 
mal vector is in radial direction, it will be a rectangle 
in the sky with angular dimensions corresponding to the 
comoving size 



CiU{z{ti) + l) X Vs^sUiziU) + l) 



(22) 



(see Eq. Q). If the angle of the normal of the string 
plane to the radial normal vector is 7^ 0, then one of the 
planar dimensions is reduced by a factor of cos{9). How- 
ever, the light travel time through the wake is increased 
by a factor of [cos(0)]~^. Thus, the wake becomes a bit 
smaller but the signal strength increases. The average 
amplitude of the polarization is given by Eq. - a 

value obtained using the average thickness of the wake. 
However, since the thickness of the wake increases lin- 
early in direction of the string motion, the amplitude of 
the polarization signal will also increase linearly in this 
direction. A second source for the linear increase in the 
amplitude is the fact that z(ti) is increasing as we go away 
from the tip of the cone (by an amount corresponding to 
one Hubble expansion time when comparing the tip of 
the cone to the end). The direction of the polarization 
vector depends on the relative orientation between the 
string and the CMB quadrupole. In Figure 2 we give a 
sketch of the signal. The amplitude of the polarization 
is proportional to the length of the arrow, the direction 
of the polarization is given by the direction of the arrow. 
Since the direction of the variation of the polarization 
strength is determined by the string and is therefore un- 
correlated with the direction of the CMB quadrupole, 
on the average the E-mode and B-mode strengths of the 
polarization signal will be the same. 

Let us now discuss the magnitude of our effect. We will 
consider the value G/i = 3 x 10~^ which is the current 
upper bound on the string tension 23] (using the as- 
sumptions about the cosmic string scaling solution made 
in these papers). Wakes produced close to the time of 
recombination inherit the ionization fraction of the uni- 
verse at that time. Taking a value of / = 10~^ (which is 
smaller than the ionization fraction until redshift z ~ 600 



Z(T )=11Q0 
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FIG. 2: The polarization signal of a single wake which is taken 
to be perpendicular to the line of siglit between us and the 
center of the string segment. The tip of the wake (position 
of the string at the time the wake is laid down) is the ver- 
tical edge of the rectangle, and the string velocity vector is 
pointing horizontally to the left The length and direction of 
the arrows indicate the magnitude and orientation of the po- 
larization vector. We have assumed that the variation of the 
quadrupole vector across the plane of the wake is negligible. 
Note that the angle between the velocity vector of the string 
and the CMB quadrupole is random. We have assumed that 
the quadrupole vector is at an angle relative to the plane of 
the wake. This determines the direction of the polarization 
vector we have drawn. The orientation we have chosen in this 
figure corresponds to an almost pure B-mode. 



[5l|), we obtain a polarization amplitude of P ^ 10^^ 
/iK which is larger than the background in the B-mode 
polarization arising from weak lensing of the primordial 



perturbations [52,|53| for 1- values of about 100. 



It is instructive to compare our polarization signal from 
a cosmic string wake with the expected noise due to the 
Gaussian fluctuations. In Fig. 3 we have superimposed 
the map of the Q-mode polarization from a cosmic string 
wake laid down during the first Hubble time after re- 
combination with a corresponding Q-mode map due to 
Gaussian noise of the concordance ACDM model. The 
string parameters are the same as mentioned in the pre- 
vious paragraph. We chose the orientation of the string 
relative to the CMB quadrupole such that the power in 
the Q-mode is half the total power. As the value of the 
CMB quadrupole we used 30/xK. To render the string sig- 
nal visible in the Q-mode map (in which the noise is much 
larger than in a B-mode map) we multiplied the string 
signal by 100. In this case the string signal is clearly vis- 
ible be eye. The brightest edge (the vertical edge on the 
right side) corresponds to the position of the string when 
it begins to generate the wake, not at the position at the 
end of the time interval being considered (the vertical 
edge on the left). Note the difference compared to the 
Kaiser-Stebbins effect in the CMB temperature maps: in 
this case the brightest edge corresponds to the location 
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IV. CONCLUSIONS AND DISCUSSION 



FIG. 3: The Q-mode polarization signal of a single wake which 
is taken to be perpendicular to the line of sight between us 
and the center of the string segment, superimposed on the 
Gaussian noise signal which is expected to dominate the total 
power spectrum. The string signal is multiplied by a factor of 
100 to render it visible by eye. The tip of the wake (position 
of the string at the time the wake is laid down) is the vertical 
bright edge of the right side, and the string velocity vector is 
pointing horizontally to the left. At the final string position 
the wake thickness vanishes and there is no polarization dis- 
continuity line. We have assumed that the variation of the 
quadrupole vector across the plane of the wake is negligible 
and that the quadrupole vector is at an angle relative to the 
plane of the wake such that the Q-mode picks up half of the 
polarization power. 

of the string when the photons are passing by it. 

Without boosting the string signal by a large factor, 
it would not be visible by eye. However, the distinctive 
lines in the map can be searched for by edge detection 
algorithms such as the Canny algorithm which was used 
to study the st ring signal in CMB temperature maps. In 
the studies of (SJ] it was found that the cosmic string 
lines in temperature maps can be picked out if the string 
signal accounts for less than 0.2% of the power. Thus, 
it should be able to easily pick out the string signal in 
the Q-mode polarization maps. In B-mode polarization 
maps the string signal would be much easier to detect. 



In this Letter we have discussed a position space signal 
of a cosmic string wake in CMB polarization maps. In 
the same way that the line discontinuities in the CMB 
temperature maps predicted by the Kaiser-Stebbins (KS) 
effect yield a promising way to constrain/detect cosmic 
strings in the CMB (see e.g. [2£-30, 32-34i], we believe 
that the signal discussed in this paper will play a similar 
role once CMB polarization maps become available. 

We have shown that a single wake will produce a rect- 
angular patch in the sky of dimensions given by equa- 
tion (22), average magnitude given by equation and 
amplitude increasing linearly in one direction across the 
patch. For a value of the string tension G/i = 3 x 10~^ 
(the current upper limit), the amplitude of the signal is 
within the range of planned polarization experiments for 
wakes produced sufficiently close to the surface of recom- 
bination. These wakes are also the most numerous ones. 
The brightest edge in the polarization map corresponds 
to the beginning location of the string, not the final loca- 
tion. Since the KS discontinuity in the CMB temperature 
map will occur along the line corresponding to the final 
position, the polarization signal discussed here provides 
a cross-check on a possible string interpretation of a KS 
signal. 

A scaling distribution of strings will yield a distribu- 
tion of patches in the sky, the most numerous ones and 
the ones with the largest polarization amplitude being 
set by wakes laid down at times close to the time of re- 
combination which are crossed by the CMB photons at 
similarly early times. 
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